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Abstract 

A practical version of the polynomial canonical formalism is developed for normal 
mesoscopic systems consisting of N independent electrons. Drastic simplification of 
calculations is attained by means of proper ordering excited states of the system. In 
consequence the exact canonical partition function can be represented as a series in 
which the first term corresponds to the ground state whereas successive groups of 
terms belong to many particle-hole excitations ( one particle-hole two particle-hole 
and so on). At small temperatures (T < 10 inter-level spacings near the Fermi level) 
the number of terms which should be taken into account is weakly dependent on 
N and remains < 10 even if ~ 10^. The elaborated method makes canonical 
calculations to be not more complicated than the grand canonical ones and is free 
from any limitations on N and T. 
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1 Introduction 



The characteristic feature of the modern stage in studying mesoscopic sys- 
tems is a possibihty to fix individual properties of investigated objects such 
as chemical composition, size geometric shape and particle number. The well 
known example of such investigations is the discovery of electronic shells in al- 
kali metal clusters [1]. The other branch of researches is connected with recent 
advances in the development of new techniques for fabricating two- or three- 
dimensional micro- and nano-structures that enables current experiments to 
investigate variations of both the exact number of electrons on such mesoscopic 
structures and their geometrical shapes with a precise control in all dimen- 
sions. Frequently experiments are produced so that investigated systems may 
exchange energy with surroundings (they arc embedded in heat reservoir), i.e. 
they are kept at an invariable temperature, however the exchange by particle 
is absent. Thus, such systems constitute a canonical ensemble, and the most 
appropriate theoretical method for their statistical description is the canonical 
formalism. 

The papers of Denton, Miihlschlegel and Scalapino [2] were the first where the 
canonical description was developed in applying to electrons in small metallic 
grains. At that time it was impossible to fix the shape and particle number of 
each grain. Therefore for constructing the canonical partition function Denton 
et al made the simplest assumption that the energy spectrum of free electrons 
in a metallic grain is equidistant ( the equal level spacing model) and two-fold 
degenerated if there is no magnetic field. Using reasonable approximations 
(M > N 1, M is number of single electron levels) Denton et al established 
an analytical expression for the canonical partition function. They showed that 
the main differences between the canonical and grand canonical values of the 
heat capacity (C) and magnetic susceptibility (x) lie at small temperatures 
T < 6 {S is the mean level spacing). Besides they found that at T 3> 5 
the grand canonical heat capacity exceeds the canonical one by /cs/2 {ks is 
the Boltzmann constant). The partition function of Denton et al takes into 
account energy variations in an applied uniform magnetic field. They manifest 
themselves in oscillations of x and C with growth of the field. Such oscillations 
are typical for charged particle systems [3] however the strong periodicity in 
the oscillations is a feature of the equal level spacing. 

It is worthwhile emphasizing the A'"-independence of thermodynamic quanti- 
ties as functions of T/S in the canonical calculations of Denton et al at M > A'", 
M,N ^ 1. It is caused by a relatively narrow layer of single particle levels 
near the Fermi energy, which mainly contributes to the partition function, and 
by the identical for any structure of the single-particle energy spectrum in 
this layer. 
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The canonical formalism was also used by Brack, Genzken and Hansen [4] 
for calculations of thermal properties of the valence electrons in alkali metal 
clusters. In this case free electron energy spectra contain considerable gaps that 
accounts for enhanced stability of clusters with the magic numbers of atoms. 
Brack et al showed that the shell structure reveals itself in the theoretical heat 
capacity. However, the difference between the canonical and grand canonical 
shell effects is quite noticeable only at temperatures which are less than or 
of the same order as the mean level spacing dp near the Fermi. Thus, both 
investigations, [2,4], point out the region of the temperatures, T < Sp, where 
the comparison of the canonical and grand canonical calculations is of the most 
interest. Besides, in the center of attention should be magic electron numbers 
(if they appear in the system) since the existence of such is a consequence 
of the appearance of an energy gap near the Fermi level that enhances the 
difference between canonical and grand canonical results. It is well known that 
T < 6f is the best regime for observing size effects in mesoscopic systems [5,6]. 
Therefore the canonical approach at these temperatures can give more precise 
theoretical information than the grand canonical calculations. 

In order to calculate the partition function Z{N, M), occupation numbers, in- 
ternal energy and other thermodynamic quantities Brack et al used recurrence 
relations connecting Z[N, M) with those for smaller N and M. Following this 
procedure they found that N{M — N) calculation steps were needed to obtain 
the final result, i.e. for M > ~ 10^ {N{M-N) ~ 10"^). Such method is still 
manageable for iV up to a few hundreds, however it could be hardly applied 
to N and M > 10^. Thus the constructing of the effectively working canon- 
ical formalism without restrictions on the particle number and applicable to 
arbitrary spacing distributions is a problem that has to be solved. 

Parallel with the exact canonical calculations in Ref. [2,4] approximate meth- 
ods have been suggested for constructing the canonical partition function of 
the normal (nonsuperconducting) systems. One of them Ref. [7] replaces the 
projection integral over the gauge variable {0 < (p < 27r) by the discrete sum 
that involves L terms. These terms depend on 27r//L (instead of ip) where the 
integer / varies from to L — 1. As shown in Sec.l this method can be viewed 
as the partial projection. It is the more precise the higher is the value of L. 
In Ref. [8] this method was employed by Frauendorf and Pashkevich to cal- 
culate sodium cluster shape and free energies. Another method suggested by 
Rossignoli consists in the evaluation of projection integral in the saddle-point 
approximation. Such approach turns out to be quite satisfactory in the case 
of not very small temperatures and particle numbers as was demonstrated by 
calculations of nuclear level densities in Ref. [9] . 

Our effective version of the canonical approach is based on the polynomial 
representation of the partition function. Last years such representation is dis- 
cussed in the literature [10,11], however our practical version of polynomial 
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calculations was elaborated many years ago for a projection method applied to 
the Bardin-Cooper-Schriffer (BCS) function for describing the particle num- 
ber conserving pairing correlations in nuclei [12]. 

There is an obvious analogy between the BCS function and grand canoni- 
cal partition function for the independent electron model {I EM). Both can 
be represented as products n(l + fs), where for the BCS function is 
(vs/tisja+a^ {vs, Us are the Bogolubov parameters, aj"a^ is the operator of the 
fermion pair creation in time conjugated states s and s) while for the canonical 
partition function fg is exp(es — A)/ kT, Sg and A being the energy of a single 
-particle state s and the chemical potential respectively. After projection onto 
a fixed particle number both functions gain the form of a iV-th order sym- 
metric polynomial [A^]. This polynomial [A^] consists of M!/A^!(M — A^)! terms 
(as before M is the total amount of employed single-particle states) and each 
term in it is the production of A" different fg, i.e. [A^] includes all distributions 
of A^ BCS pairs or A^ electrons over M states. For large M and A^ the total 
number of terms in A^ may be enormous, however a hierarchy of terms can be 
established according to their surviving at G — > (G is the pairing strength) 
or T" ^ 0. The first term of this hierarchy corresponds to the independent 
fermion ground state in which all single-particle states are occupied up to s^. 
The second group of terms corresponds to simplest excitations: in the BCS 
theory it is a particle pair a^a^ above ep and a hole pair 0^05 beneath ep^ for 
free electrons in a cluster it is one particle-hole excitations. The third group of 
terms corresponds to a shift of two BCS particle pairs or two cluster electrons 
above Sp. This alignment of many particle- hole excitations continues till all 
pairs or cluster electrons are lifted above Sp- 

The common property of the projected BCS and the canonical partition func- 
tions is mutual independence of probabilities of particle and hole excitations. 
This probability depends only on the energy of the single-particle level on 
which a BCS fermion pair or a cluster electron is located. For the BCS case 
it is (vg/us) for levels with Eg > ep while for levels with < e^? it is {ug/vg) , 
for cluster electrons it is exp(A — Eg/kT if Eg > Ep, or explEg — X)/kT if 
^ ^F- Thus summing over particle and hole excitations can be performed 
independently. The contribution of each group of n-particle excitations to [A"] 
includes a sum running distributions of n electrons (here we mean cluster elec- 
trons) over particle levels (with energies above Ep). These sums can be also 
represented as symmetric polynomials, now they are of n-th order (n < A^). 
Analogous polynomials arise for hole excitations, but in this case electrons are 
distributed over hole levels (below and including ep). Though the convergence 
of appearing series is absolute its rate is determined by T/5p. Practically for 
T/Sp < 1 they include not more than 5 components (i.e. up to 5 particle- 
5 hole excitations). Our calculations for the pairing problem have indicated 
high effectiveness of the polynomial representation of the A^-projected BCS- 
function [12] and this paper is devoted to the adaptation of the polynomial 
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representation to the canonical partition function of normal mesoscopic sys- 
tems. 

The material is arranged as following. In Sec. 2 separation of the full configu- 
ration space into two parts, inside and outside the Fermi sphere, is realized to 
construct the canonical partition function Zjv as a symmetrical polynomial of 
iV-th order. Z]\j is expanded in polynomials of lower orders. These series are 
composed of products of "particle" and "hole" polynomials, the order of which 
points out how many particles and holes take part in excitations allowed for 
in Zjv. In Sec. 3 these polynomials are employed to one- and two-body den- 
sity matrices and on this base the canonical expressions for the heat capacity 
and magnetic susceptibility are compared to the grand canonical ones. Sec. 4 
collects recurrent relations for polynomials wc apply in calculations. Sec. 5 
shows that in a wide range of temperatures and particle numbers the elabo- 
rated method possesses high convergence and can be used for various types 
of single-electron energy spectra. The conclusion is given in Sec. 6. Algebraic 
details for the equal level spacing model and the exact canonical partition 
function of the model are considered in Appendix. 



2 "Particle" and "hole" symmetric polynomials and the canonical 
partition function 

If a system is composed of a fixed number of constituent particles and 
possesses a discrete excitation spectrum the thermodynamic properties of the 
canonical ensemble iCE) of such systems can be described by using the par- 
tition function Zjv 

Zm = Y. exp{-f3E^{N)}, f3 = T'' (1) 

a 



The index a labels energy states of this system Ea{N), a = concerns the 
ground state. The calculation of Ea{N) for many body systems is a compli- 
cated problem and in this paper we limit ourselves by the I EM which in many 
cases is a good start approximation. Then, assuming that a set of single elec- 
tron energy levels is established for a static mean field governing the electron 
movement, the values of E^ are readily found in I EM. 

The ground state energy Ea{N) is the sum of A^ single- electron energies up to 
the Fermi energy (ep) and the excited state energies are obtained by addition 
to Eq of several single-electron particle energies Sp > ep) and by subtraction 
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of the same amount (as N is fixed) of hole energies {sh < Sp) 

Eo{N) = j:{es-X), (2) 

s=0 



E^{N) = Epip2 ...pk, hih2 ...hk^ Eo{N) + + Sp^ + . . . + Sp^ - 

-£hi - - ■ ■ ■ - ^hk- (3) 

We shall treat A in Eq. (2) as a chemical potential, however, in contrast to 
the grand canonical ensemble (GCE) where A is calculated to allow for the 
particle number conservation on the average, the choice of A for as a 
certain level for counting energies is merely a matter of convenience (in our 
calculations A is taken as the grand canonical chemical potential at given T). 
Substituting Eqs. (2,3) into Eq. (1) one can see that the canonical partition 
function Zjv is a symmetric polynomial in Qs 

Zn^[N]^ E f[{qs)i, qs^exp[-P{e,-X)]. (4) 

{distr.} i=l 

The symbol {{distr.} in Eq (4) implies that all A^-electron distributions over 
all single electron levels have to be taken into account. 

Another way to introduce Zn as a symmetrical polynomial consists in em- 
ploying the statistical operator Dn projected onto a fixed particle number 
N: 



Dn = -^Pn n exp{-p{es - A)a+aJ, (5) 

s 



1 

Pjv = — I dip exp{iip{N - TV)}, ^ = E 



af, tts are the creation and annihilation electron operators. The normalization 
of Dn, Zn, is calculated with unrestricted electron basis {Fock space) \) 

Zn^^I dipe^^"" Z{^y, (6) 



Z{^) = (I Y[exp{[-I3{es - A) - iv^]a+a,} |) = 

s 

^l[{l + exp[-P{e,-X)-i^]}. (7) 
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The Fourier series for in powers of e includes as amplitudes symmetrical 
polynomials [n] of the same type as in Eq. (4) 



ZH = l + ^[n]e— '^j (8) 

n=l 

[i] = E?- [2] = E?^?*; E n(a (9) 

s s<t {distr.} i=l 

The projection in Eq. (6) isolates the polynomial with n — N that again leads 
to Eq. (4). 

If integrating in Eq. (6) {Z{ip) from Eq. (8)) is replaced by summing over 
discrete variable 27r//L as was suggested in Ref. [7] 



'-'^ 

= ^ [l _ e-^-] {l - e^2./(Ar-n)/L-e)| 1^^^^ 

1 I AT-n h/cL; A; = 0, 1, . . . 
= < ' ' ' ' ' ' (10) 

0\N-n\^kL + u; I < u < L - 1, 

then the number projection does not eliminate completely components in the 
grand canonical Z which correspond to n 7^ A^. However the more is L the 
larger is the difference between A*" and such n's (not less than L). 

The expectation values of the one- or two body density matrices can be also 
expressed by means of polynomials 



Us = ajas, 



s 

Us = (I UsD^ I) = ^)^^^i (11) 

nst = (I UsUtDN I) = gsg J^^^^^^'* , s^t. (12) 

If state s is degenerated and can be occupied by several electrons, in polyno- 
mials there appear powers of the corresponding amplitude Qg. For such cases 
we introduce Ugs' 

Uss^il flsfltDN \)\et=es = ^5^^~~pvf^- (1^) 
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Henceforth the notation [iV]o,5 imphes the exception of single electron states 
a,b . . . out of the full set of states over which n particles are distributed. 

For comparison we write out here the GCE results for the same quantities: 
Z{GCE)^ll{l + q,) = l + Y.[n]. (14) 

s n=l 

The maximum n in Eq. (14) is determined by the amount of the single electron 
states under consideration. 

{\hsD{GCE) |) = /, = l/(l + g;i), (15) 
{\nsntd{GCE) |) = s ^ t, (16) 

fs is the Fermi occupation numbers. 

Eq. (12) shows that for the canonical ensemble the two-body density matrix 
Ust is not reduced to a product of two one-body ones Usfit as it takes place in 
GCE, Eq. (16). Thereby the canonical description allows for correlations in 
electron movement in spite of the starting assumption of the I EM. 

The method we are about to develop here for computations of polynomials 
entering into Eqs. (4), (11)-(13) takes advantage of the general property of 
symmetrical polynomials: each polynomial given in a space of variables 
can be represented as a decomposition of polynomials of the same and lower 
orders defined in subspaces of smaller dimensions. That can be viewed as a 
consequence of the Clebsch-Gordan decomposition apphed to the basic vectors 
of the symmetric group irreducible representations. If the original space C 
spans a set of single electron states which can be divided into two subsets 
{C = A® B) m accordance with energy or any other quantum numbers then 
an initial polynomial [n] ^'"^ can be represented as following 

n 

[nf) = ^ [n - [vf^ ■ [0](^) = [0](^) = 1. (17) 

Superscripts in this equation denote the spaces in which polynomials are de- 
fined. The order of each polynomial [n — z/]'-^-' or [u]^^^ in Eq. (17) cannot 
exceed the dimension of space A ox B respectively. 

The first realization of this decomposition is dividing the whole single elec- 
tron energy space into two parts: below and above the Fermi level {F) that 
drastically simplifies the calculation of the polynomials and hence, as it will 
be shown in Sec. 5, enables the canonical approach to be employed in a wide 
range of temperatures and particle numbers. 
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To implement such dividing we have to exphcitly allow for possible degeneracy 
of electron states and in particular the degeneracy of the Fermi level, dp- If 
dp > I the occupation rate of F dictates three ways to factorize of [N] 



Zm = [N] = Z^'^^Zm (18) 



The first case corresponds to the completely filled Fermi level at T = {np — 
dp, Up is the particle number on F-level). Then the first factor on the right 
side of Eq. (18) is the polynomial [iy]*-^^ where space A of the dimension N 
involves all states below and including F. Therefore this polynomial consists of 
only one term and is the product of all Qg with s running states < F. Factoring 
out this term in Eq. (18) leads to a series in symmetrical polynomials: 



N 

= n = [-pEom , = E N (19) 

s<F n=0 

IoI=[o] = i, 

dg is the degeneracy of state s {dg is equal to the maximum electron number 
on this level); [n] and [n] are n-th order polynomials, the former is a "particle" 
polynomial defined in variables qg = exp[(3{X — Sg)] = fg/{l — fg), Sg > Sp 
while the latter is a "hole" polynomial in q^^ — exp [/3{£t — A)] = (1 — ft)/ ft: 
^t < fs is determined in Eq. (15). Thus, A being chosen between F and 
F + 1 (the level just above F), both [n] and [n] are functions of quantities 
which are < 1. This choice of A and exponential energy damping of variables 
Qg and q^^ provide the prompt convergence of the series in Eq. (19) especially 
at small temperatures. 

The second case is np < d//2 at T = 0. Now the factors Z^^^ and Zjv are 

F-l _ N-np 

= n {<lsf-. ^iv = E ¥][nF + n], (20) 

s n=0 



F — 1 is the topmost state just below F. 

The third case is dp/2 < np < dp a,t T — 0, then 

F _ ^ 

^^°^=n(5*)''; Z^ = J2idF-np + n][n]. (21) 

t n=0 



9 



3 Canonical magnetic susceptibility and heat capacity 



The polynomials [n], [n] and those with omitted states determine one- and 
two-body density matrices (Eqs. (11-13)). Their explicit form for the rip — dp 
case is following ( is defined by Eq. (19)): 



1 - (grV^iv) EtiN[n - 1],, iis<F ^^^^ 
(9./^iv)En=i[n-l]sR, if s>F; 

Ust = si^t] (23) 

Qt - Qs 



+(g.-V^iv)Et2W[ri-2L, if.<F 

(g.7^iv)Et2[^-2]..R, ifs>F. 



(24) 



These density matrices and analogous ones that correspond to Eqs. (20), (21) 
determine Xcan and Ccan: 



Acan Acan ~ Acan ) 



(25) 



Superscripts P and D mean the paramagnetic and diamagnetic parts of %. 



Acan 



V 



0£ 

+ E^^.(4-l)(^)^(n..-n^) 



du; 

-F2^4(it(-^) {-^){nst - nsUt) 



(D) _ 



(26) 
(27) 



C'can/ ks — P'^ 



s 

+ ds{ds - 1)(£. - A)2(n,, - nl) 

s 

+ 2Ydsdt{es - A)(£t - X){nst - nsUt) 



s>t 



(28) 
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The correlations caused by the particle number conservation formally reveal 
themselves in the appearance two in Eqs. (26,28) terms containing two-body 
density matrices (n^j and Ugs)- Besides, due to the correlations, i.e. due to the 
dependence of the occupation probability of a level on the occupation of the 
other levels, the Fermi occupation numbers fs differ from the canonical ones 
Us- The grand canonical Xgrand and Cgrand can be obtained from Eqs. (26,28) 
by replacing by fs and removing terms proportional to {ugt — nsUt) and 
{nss—n1). Both these factors now vanish since in the absence of the correlations 
in the grand canonical ensemble two-body matrices are factorized: {ust = ngrit) 
and {uss ^nl). 



Xgrand Xgrand Xgrand' 

(29) 

xZ.-^PEds(^^Jfs{l-fs), (30) 



(D) _ /o-|N 
Xgrand ~ T/ * fl, ,2 ' 



C^.nJkB = P^Y.ds{es - A)/.(l - fs){e, - A - /3— ), (32) 



d\ _ Eds{es-X)fs{l-fs) 
dp PEdsfsil-fs) 



(33) 



Eq. (32) for Cgrand is written at the condition that the average value of N 
is conserved. In this case the chemical potential is subject to temperature 
variations, term in Eq. (32) pdX/dp. This term is absent if measurements are 
performed at a constant A. An analogous term in the susceptibility is dropped, 
since it is vanishingly small. 



4 Calculations of the polynomials 



Polynomials entering into Z^, Xcw and Ccan can be expressed through Qs and 
q^^ with help of recurrent procedure based on Eq. (17). If in this equation the 
partition C = A Q) B is such that space B consists of only one state s then 
[u]^ = ^ly.iQ's + Su.o and polynomial [n — v^^^ in space A is in fact \a — v\s 
since it is defined in the whole space C with state s omitted. Eq. (17) can be 
therefore rewritten as 

[n] = [n - \\sqs + [n],. (34) 
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The sum rule gives the second equation for successive calculations of polyno- 
mials: 

-Y^[n-l]sdsqs = [n]. (35) 



The meaning of Eq. (35) is illustrated by its apphcation to Eq.(lO) 

J2[N - lldsQs = [N]{\ NDn I) = [N] ■ N. (36) 

s 

Here s enumerate all single electron states. 

The lowest order polynomials being found straightforwardly 

[0] = [0]i = l; [l]^Y.dsqs, (37) 



higher order ones can be determined by using Eq. (34), that gives the well- 
known Newton identity (35) 

[n]^-f:[n-u]a^{-ir-\ (38) 

a^^^d^Qs^ (39) 



Eq. (38) can be used to express [n] only through a,^, Eq. (39), such formula is 
given in textbooks on algebra, but it is not so convenient for calculations as 
Eq. (38). 

n 

M. = E[^-%r(-ir- (40) 

i/=0 



If Qs 7^ Qt for states s and t then 



V 



= (?tNt - 9sNs)/ {it - Is) (41) 

^{[n + l]t~[n+l],)/{qs-qt). (42) 
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If state s is degenerated {ds > 1) then it is possible to define [n\ss, Eq. (13), 
as following 

N.. = E[^-H(-ir?:(^+i) (43) 

i/=0 

or by using recurrent relations 

[0],, = 1, [nls = [n\s - qs[n - 1],, (44) 

Eqs. (38), (40-43) are sums of sign-changing terms. Instead of them one can 
use sums with positive terms 

N^''^^ E [^llf^^■[^l2f^^■...■[^lnf-\ (45) 

that can be obtained from Eq. (17). Superscripts in Eq. (45) assume that 
C — Si® 82® . . . (B Sn where each subspace Si embraces d{Si) single electron 
states, therefore /Xj < d{Si). If these states are degenerated, any polynomial 
in Eq. (45) is calculated analytically 

= C'^^'^expi-P^iiss - A)}, ) = "^^^^ j . (46) 



is a binomial coefficient. Eqs. (45), (46) are appropriate to calculate the 
derivatives of [n] with respect to Qs 

M^d,[n-lU ^ ^ ds{ds - l)[n - 2U. (47) 



If each subspace Si is nondcgcncrate, i.e. each yUj = or 1, Eq. (45) is trans- 
formed into the sum over distributions similar to Eq. (4). 

When a magnetic field is applied to a system in which spin and orbital mo- 
menta are decoupled each electron spin is aligned along or against the field B. 
That gives rise to the Zeeman splitting of the single electron levels 

e,^e,{B)±u;, u;^fj,*s^B, (48) 



IJ,*^ is the effective Bohr magneton and g is the effective gyromagnetic ra- 
tio. es{B) incorporates magnetic effects on the spatial electron density. The 
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spin splitting can be explicitly separated in Zn by means of dividing the 
whole space C of single electron states into two subspaces: C = (7*^+-* © C^^\ 
spans states with positive magnetic spin contributions H-a; to the energy, 
Eq. (45), C*^~) corresponds to —u. Such partition causes the decomposition of 
each polynomial: 

[n]^Y.[n-v\M[v\M- v^^mm{n,G). (49) 

v=0 



G is the whole amount of single particle levels which were twice degenerate 
in the absence of the field. Polynomials [t^]*-^-* include factors expl^Puju] that 
allows one to reduce Eq. (49) to the form 



[n] — exp{— ^ exp{2Pu;u}[[n — 



(n-5n)/2 

Y: [2 - {1 - Cn)Su.o] 
1^=0 



i/=0 



+ V 



— V 



X 



cosh \fiuj(2v + , 

en = [l-(-)1/2. 



(50) 



The polynomials [[(n±^„)/2±z/]] in Eq. (50) are composed of independent 
of the spin magnetism: = exp{—l3[es{B) — A]}, es{B) is given by Eq. (48). 
Eq. (50) is employed in Appendix to find Zn for the equal level spacing model. 
Eqs. (35 - 50) are written out for "particle" polynomials [n] for which space G 
stretches from level (F+l) up to the topmost level M. The same equations are 
applied to "hole" polynomials [n] provided variables qg with {es{B) — A) > 
are replaced by q^^^ with {st{B) — A) < 0. The level space for [n] ranges along 
all single-particle space beneath and including level F. 



5 The potentialities of the polynomial method 



As mentioned in Introduction the best temperature regime to reveal the dif- 
ference between canonical and grand canonical description is T < Sp- At this 
condition the convergence of series over n for Zj^, Eqs. (19)-(21), is attained by 
taking into account as low as a few terms {n <^ N). However, in this section we 
will show that the polynomial method can be applied in a wide temperature 
range since these series possess the absolute convergence. Moreover, even for 
such high temperature as T ~ 506f real values of n do not surpass ^max ~ 100 
for iV ~ 10^ i.e. TT-max far smaller than iV, and the most part of series for 
Zn, which in the general case stretch to N, does not contribute to Zj^. 
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Firstly we shall accomplish estimates of the convergence analytically and ex- 
ploit for this purpose the equal level spacing model that is reasonable approx- 
imation to a real system at high temperatures T > dp when details of possible 
shell structure are averaged out. To avoid the double degeneration of levels 
we dispose the system in a weak magnetic field B to actualize the Zeeman 
splitting, Eq. (48). If a; = A/4 {qs = 2) the levels arc cquidistantly placed as 
before with the level spacing A/2. Here and in Appendix A is used for the 
level spacing at -B = 0. In this case the polynomials entering into 

ZN = Y^z{n), z{n) = [n][n], (51) 

n=0 



are defined in a space of "particle" (qs) and "hole" {q^ ^) variables, Eq. (4), 
reduced to powers of q (here we put X — ep)'- 

g = exp(-l/t), t = 2T/A (52) 
qs ^ exp[-P{eF+s - Sf] ^ q^ (53) 
q^^ = exp [-p{eF - SF-t)] = (54) 

As shown in Appendix polynomials in such variables can be found analytically, 
Eq.(A3). li M > N ^ 1 these expressions are simplified 

H = g'^("-i)/2 f[{l- g^)-S [n] = g^p. (55) 
t/=i 



At small T when t -C 1 the series for Zjv is exhausted by few first terms since 
g -C 1. However the convergence falls if T increases giving rise to g ~ 1. At 
such temperatures the polynomials [n] are increased with n up to some no 
after which they begin decreasing. An estimation of no is obtained from the 
condition: 

'-^^^<1, ^%±il<l, (56) 
z{no) z[no) 



Substituting Eq. (55) into Eq. (56) gives 

no ~ Un(l + exp{l/20) ^ 0.7t. (57) 



Thus no increases with t and practically does not depend on N ii no <^ N. 
The upper limit of summation, n^iax, can be found from the condition 

^(^max ~t~ 1) ^ Q 
-^(^max) 
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This equation yields nmax ~ 2no that is demonstrated by Fig. 1 where the 
results of the exact calculations of Z^r for the equal level spacing model are 
given. It displays the ratios of Zjv calculated with few terms of the full sum 
in Eq. (51) to the exact value of Z^. At a fixed temperature the bend points 
for N — 10^ 10^ practically coincide. These points correspond to no and the 
saturation sets in even earlier than at 2no. 

As shown in Fig. 2 the values of hq at a fixed reduced temperature ( T/5f) 
depends on the type of the electron confinement. The energy single-electron 
spectra of the equal level spacing model and that of the anysotropic oscillator 
are similar. Therefore similar is the convergence of the corresponding Zpf- 
series. For a spherical cavity the spectrum reveals shell gaps. In consequence 
the convergence turns out to be somewhat higher. 

At high temperatures f/ ~ 1 — 1 ft and 1 — ~ n/t, therefore the n-dcpcndence 
of any term z{n) in Eq. (51) at n ^ 1 can be obtained from Eq. (55) 



that proves the absolute convergence of Zjv-series. 

Asymptotic values of Zjv at high temperatures can be found if summing over 
single-particle states is replaced by integrating with that part of the level 

density p{e) which is monotonously varying with energy e. At such T the 
fluctuation part of p{e) (shell correction) practically brings no contribution. 
Along with this approximation two other conditions are employed: T <^ Sp 
and > 1. 

in I EM defined as an integral, Eqs.(6),(7), can be represented as following 



{t) 



(58) 





(59) 







A= ^ln{l + exp[/3(£, - A) +i<^]} 



s<F 



A 




(60) 







B = ^ In {1 + exp[ - /?(£t - A) + ]} 



OF 



oo 




(61) 
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For N ^ 1 the upper limit A in Eq. (60) is extended to oo. Then taking 
into account only first term (independent of (3) in the expansion of functions 
p(A ± yP~^) one arrives at 



A + B = p{X)l3-^ y" In (l + 2 cos ipe-^ + e'^y) dy 



Z (tT -(^ ), < TT, 



Thus Zjv gains the form 



(62) 



= exp {-l3{Eo - \N) + ^p{\)/Qf3] 
e = -/d^exp[-p(A)/3-VV2] 



(63) 



As p{X)f3 ^ is supposed to be » 1 the upper limit in the integral in Eq. (63) 
can be taken as oo: 



1/2 



.27rp(A) J 



exp 



-p{Eo - XN) + 7rV(A)/6/? 



(64) 



The value of Zg^and is found from Eq. (63) if ^ = 1 
^grand = exp \-P{Eo - XN) + ^p{X)/&(3 . 



(65) 



In Eq. (65) A is regarded as a function of T defined by the particle number 
conservation. 



Now one can see the difference between Ccan and Cgrand 



,d'^\nZN 



4Tp(A)7rV3 - kB/2; 



(66) 



C'grand — ^s/? 



,dU 



grand — '"at-' : 

d In -Cgrand i ^-1 ^ 1^ -Cgrand 

\x=const — , 

Cgrand = k%Tp{X)n^/3. 



(67) 
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In Eq. (67) the condition T/sp <^ 1 and the relation 

p{X)dX/d(3 ~ [5-^{dp/d\)Ti^/2, 

are employed. Thus the difference between Ccan and Cgrand is the same as it 
was found in Ref. [2] for the equal level spacing model. Eqs. (66,67) show that 
the same difference exists for any level density. 

As shown in Appendix the exact canonical partition function for the equal level 

spacing model turns out to be much more complicated than the expression 
given in Ref. [2] in approximation M > ^ 1. However the difference reveals 
itself mainly at small particle number and T/A > 1. This is demonstrated in 
Figs. (3), (4) where Ccan is displayed v.s. A^ and T. 



6 Conclusion 

We have shown that the canonical partition function Zj^ which in general case 
is a symmetric polynomial of the A^-th order (in variables of exp{/9(A — e^)}) 
can be expanded in polynomials the orders of which range from to A^. Such 
representation of takes into account that in the independent particle model 
the full configuration space is naturally divided by the Fermi surface and each 
excitation in this model is generated by particle- hole pairs the amount of which 
range from to A^ as well. Therefore the polynomial expansion of Zjv is a sum 
over products of particle and hole symmetric polynomials the order of which 
vary in the same manner as the amount of pairs in the excitation. 

In consequence the canonical averages such as the occupation numbers, heat 
capacity and magnetic susceptibility arc expressed via similar sums involving 
the same polynomials and others defined in narrower spaces (where a part of 
states is omitted). We have suggested recurrent relations reduced the calcula- 
tions of n-th order particle (hole) polynomials, corresponding to distributions 
of n particles (holes) over states above (below) the Fermi level, to first order 
polynomials. 

Though the sum Z^v stretches to A^-th order polynomials, its convergence is 
so high that this version of the exact canonical approach can be employed for 
mesoscopic systems in a wide range of temperatures and particle numbers. 
Moreover , as only few first terms bring in a real contribution to canonical 
averages in the temperature range of several level spacings near the Fermi 
level, the elaborated method is especially effective for studying structure and 
size effects in the temperature regime which is most appropriate to this goal. 
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A The canonical peirtition function of the equal spacing model. 



In this model each single electron level in the applied magnetic field B gains 
the Zeeman splitting 

ef = sA±u (A.l) 



A is the level spacing, s = 0,1, . . . , M — 1, where M is the amount of the 
electron levels, M > N. N is the electron number. Thus each polynomial in 
Eq. (50 is defined in the space of Qs (here we set A = 0) 

q, = q', q = exp{-(3A), (A.2) 
that allows a polynomial [[k]] to be calculated analytically 

--q''^''-'^/^Q{M,k), (A.3) 



k 

Q{M, k)^5k.o + (1 - M n(l - g^+'-'^V(l - Q^) 

K=l 

Eq. (A.3) is straightforwardly found for all k at small M and then extended 
for any M by induction. Substituting Eq. (A.3) into Eq. (50), where n has to 
be replaced by A^, one obtains 

= [N] = exp{-(3Eo{N)}Z^, (A.4) 



2 



u=0 



xQ(M, + u)Q{M, - v) cosh{/3(2i. + ^iv)}, 

1 -(-!)» 
= ■ 



Eq{N) is the ground state energy in this model T — B — 0; Eq{N) — 
A[Ar(7V-2)+eiv] /4. 

Eq. (A.4) is the exact partition function. In order to crossover to the expression 
given in Ref.[2] two approximations have to be done. Firstly the amount of 
the levels M has to be infinitely large, M — > oo, then 

k 

Q{M, k) ^ g(oo, k) = Sk.o + (1 - Sk.o) n (1 - qT'- 

K=l 
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Secondly the values of g = exp{—(3A) must be small enough whereas the 
particle number ^ 1 then the dependence Q on z/ in Eq. (A. 4) can be 
dropped 



References 



[1] W. D. Knight, W. A. Heer, and K. Clemenger Solid State Commun. 53 445 
(1985) ;W. A. Heer, W. D. Knight, W. D. Chou, and M. L. Cohen Solid State 
Phys. 40 93 (1987). 

[2] R. Denton, B. Miihlschlegel, O.J. Scalpino, Phys. Rev. Lett. 26, 707 (1971); 
Phys. Rev. B7, 3589 (1973). 

[3] N. K. Kuzmenko and V. M. Mikhajlov, Phys. Lett. A296, 49 (2002). 

[4] M. Brack, O. Genzken and K. Hansen, Z. Phys. D21, 65 (1991). 

[5] W. P. Halperin Rev. Mod. Phys. 58 533 (1986). 

[6] J. A. A. J. Perenboom, P. Wyder and F. Meier, Phys. Reps. 78, 173 (1981). 

[7] G. H. Lang, C. W. Johnson, S. E. Koonin, and W. E. Ormand, Phys. Rev. C48, 
1518 (1993). 

[8] S. Frauendorf, V. V. Pashkevich, Surface Rev. and Lett. 3, 241 (1996). 
[9] R. Rossignoh, Phys. Rev. C51, 1772 (1995). 

[10] H.-J. Schmidt and J. Schnack, arXiv: |cond-math/0104293| 2001; arXiv: 
icond-math/0209397, 2002. 



[11] F. Philippe, J. Arnaud, L. Chusseau, arXiv: |math-ph/0211029| 2002. 

[12] N.K. Kuzmenko, V. M. Mikhajlov, Buh. Acad. Sci. USSR, Ser. phys., 37, 1911 
(1973); Sov. J. Part. & Nucl., 20, 830 (1989); Phys. Rev., C51, 2826 (1995). 



20 




21 




Fig. 2. Influence of the temperature on the convergence of series Zjv for different 
model spaces. Zj<;{k) is defined in the capture to Fig.l. 5 is the mean level spacing 

near the Fermi level. The right panel shows fragments the single-electron spec- 
tra near the Fermi level for these model spaces: (a) the equal level spacing; (b) an 
anysotropic oscillator potential with frequencies oOx/'jJy = 1-33, uJx/ijJz = 1-55; (c) a 
spherical cavity with R = vqN^^^. 
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Fig. 3. The canonical heat capacity vs temperature. The bold curve represents the 
results of Ref. [2] and coincide with the exact canonical results if A'^ > 10^. For 
smaller N the high temperature (T > 1.5(5) behavior of Ccan (the curve for N = 20, 
M > ION, AI is the number of spin degenerate single-electron levels) differs from 
predictions of Rcf. [2]. The curve for N = 20, N = M shows that high temperature 
behavior of C in small electron ensembles depends on M if M < lOA'^. 
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Fig. 4. The canonical electronic heat capacity vs the particle number in the equal 
level spacing model {b is the level spacing) . Solid lines represent the results obtained 
by our method. The dashed line represents the result of Ref. [2]. At T = (5 the exact 
canonical results coincide with the results of Ref. [2] . 
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